Emergence of central mode in the paraelectric phase of ferroelectric perovskites 
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THz-range dielectric spectroscopy and first-principle-based effective-Hamiltonian molecular dy- 
namics simulations were employed to elucidate the dielectric response in the paraelectric phase of 
(Ba,Sr)Ti03 solid solutions. Analysis of the resulting dielectric spectra suggests the existence of a 
crossover between two different regimes: a higher-temperature regime governed by the soft mode 
only versus a lower-temperature regime exhibiting a coupled soft mode/central mode dynamics. 
Interestingly, a single phenomenological coupling model can be used to adjust the THz dielectric 
response in the entire range of the paraelectric phase (i.e., even at high temperature). We con- 
clude that the central peak is associated with thermally activated processes, and that it cannot be 
discerned anymore in the dielectric spectra when the rate of these thermally activated processes 
exceeds certain characteristic frequency of the system. 
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It is well known that the static permittivity of fer- 
roelectric materials is related to frequencies of all polar 
phonon modes through the Lyddane-Sachs- Teller formula 
|l|. Near the phase transition, however, an additional 
low-frequency mode has to be often taken into account 
- the so-called central mode A generic rea- 

son for this additional Debye-type excitation seems to 
be large-amplitude fluctuations between quasi-stable off- 
center ionic positions. Existence of such intrinsic central 
mode could be very clearly demonstrated, e.g., below the 
cubic-tetragonal phase transition Tq of BaTi03 Q • 

Similar central mode (CM) is also known to exist in 
the paraelectric phase. Phenomenological theories of the 
paraelectric CM have been developed by several authors 
[246|- All these approaches lead to a coupled relaxator- 
oscillator dielectric response. However, an important 
question has been left open so far: whether the CM per- 
sists up to the highest temperatures, or rather it pro- 
gressively disappears, or whether it disappears at some 
well-defined temperature Tqm(> Tq). 

Unfortunately, it is much more difficult to obtain a 
clear-cut experimental evidence for the dielectric CM in 
the cubic perovskite phase @, 0, @, @] • The characteris- 
tic frequencies of the soft phonon-oscillator and CM in 
KNb03 and BaTi03 are so broad and close together that 
they can hardly be disentangled. The difficulty of the 
experimental analysis of the dielectric spectra of KNb03 
and BaTi03 is caused, at least partly, by intrinsically 
large phonon damping factors, related to the fact that 
T c is quite high in these compounds (about 700 and 
400 K, respectively). Therefore, it is of interest to study 
the CM also in the ferroelectric compounds in which the 
phase transition occurs at lower temperatures, which can 
be achieved for example by mixing BaTiC>3 and KNb03 
with suitable incipient ferroelectrics, namely SrTi03 and 



KTa03, respectively. 

Here, we describe a combined experimental and 
theoretical study of the technologically-relevant mixed 
Ba^Sri-^TiOs system (BST) aimed to establish the char- 
acteristic temperature trends in the model parameters 
of the relaxator-oscillator dielectric response in the po- 
lar perovskites. It turns out that this phenomenological 
theory allows to "easily" understand the existence of the 
temperature Tcm, at which the CM in the dielectric spec- 
tra is appearing or disappearing. 

Basic dielectric properties of BST solid solutions and 
the concentration diagram of this system were thoroughly 



discussed, see e.g. Refs.llOl. Illl and references therein. 
Experiments discussed below were carried out with a set 
of high-density BST ceramics with Ba:Sr ratio ranging 
from x=0 to x=I, prepared by methods described else- 



where (see Refs. [12l . ll3l |). Low-frequency permittivity ob- 
tained from standard dielectric measurements (10 kHz) 
in the paraelectric phase was fitted to a Curie- Weiss law 



eiO)- 1 = (T-T )/C 



(1) 



As expected 0, [ijj, Curie constants were of the order 
of 10 5 K for all concentrations, while the extrapolation 
temperature Tq(x) turned out to be strongly dependent 
on the concentration (see the inset of Fig. la). 

These findings may suggest that at a fixed tempera- 
ture interval above the extrapolation temperature To, the 
static permittivity is roughly independent of x. However, 
this is not the case for the THz-range dielectric spec- 
tra of BST ceramics. In Fig. la, we compare dielectric 
loss spectra obtained similarly as in Refs.[l2|, [l_3| from si- 
multaneous fits to time-domain THz transmission and 
far-IR reflectivity spectra. Although all these spectra 
are obtained at 100 K above the corresponding tempera- 
ture Tq, we can see a clear crossover between two rather 
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FIG. 1: (Color online) Dielectric loss spectra of BST ceramics 
determined (a) by combined IR reflectivity and time-domain 
THz transmission spectroscopic techniques and (b) by MD 
simulations using the effective Hamiltonian of Ref. [111 ]. The 
spectra are obtained at 100 K above the extrapolation tem- 
perature To of the Curie- Weiss law, indicated in the insets 
by full line. Point symbols in insets indicates temperature- 
concentration coordinates associated with the spectra shown 
in the main panel. 



distinct regimes: a broad-band regime for BaTiCVrich 
compounds, and a narrow-band regime for SrTi03-rich 
compounds. The narrow response band can be easily ad- 
justed with a single damped harmonic oscillator (DHO) 
model, while the broad band suggests a more complex 
spectral shape, such as that of the coupled relaxator- 
oscillator model. 

In order to prove the intrinsic origin of this spectral 
broadening phenomenon, we have also calculated [151417 
dielectric spectra of disordered BST solutions by con- 
ducting molecular dynamics (MD) simulations using the 
effective Hamiltonians of Ref. |Tl|- as previously reported 
in Refs.0,B Chemical disorder in the Ba and Sr ion sub- 
lattice has been included directly by specifying randomly 
chosen Ba and Sr ion positions within a simulation super- 
cell representing a 12x12x12 perovskite formula units of 
BST, assuming periodic boundary conditions. To verify 
that the results are independent of the particular con- 
figuration, simulations were also carried out for several 
different configurations. 

The resulting spectra (see Fig. lb) are quite 
monotonously varying with the average Ba concen- 
tration in the simulated supercell, even though the Ba 
ion distribution is obviously different for each concen- 



o 

X 



-CM 





/ 


simulations: 


j\ 


270 K _ 




700 K 


\ SM 


fit: 


— T=270 K 


V(\ 


T=340 K 




T=700 K 


- C M , ' \\ \ 













50 100 150 200 250 300 

Wavenumber (cm' 1 ) 

FIG. 2: (Color online) Dielectric loss spectra of 
Bao.sSro.sTiOa determined by MD simulations corresponding 
to T = 270 K (at 25 K above the theoretical extrapolation 
temperature To), at T = 345 K and at T — 700 K. Smooth 
continuous lines are fit to the model of coupled oscillator- 
relaxator model of Eq.[2] (fitting parameters are shown in 
Fig. 3.) 



tration. This also confirms that the observed trends 
are driven by the overall concentration, rather than by 
the degree of the occupational disorder. As consistent 
with the experiments, the predicted static dielectric per- 
mittivity above Tc obeyed the Curie- Weiss law (Eq. 1). 
Moreover, simulated loss spectra for temperatures about 
100 K above Tq indeed confirmed a similar broadening 
for BaTi03-rich compounds as in the measurements (see 
Fig. lb). 

Despite of this considerable broadening, spectra of 
Fig. 1 do not show a very clear 2-maximum dielectric loss 
response, such as the one observed in the Ai spectra of 
the tetragonal phase of BaTiOa single crystal [7J . We ex- 
pect that CM and soft mode (SM) components are much 
better separated in the spectra within the first few tens of 
K above Tq- However, in this temperature range the THz 
measurements were not feasible even with very thin ce- 
ramic samples because they become too opaque for THz 
radiation. 

Motivated by this difficulty, we have calculated dielec- 
tric spectra of Bao.sSro^TiOa (BST50) in a broad range 
of frequencies and temperatures down to the vicinity 
of the ferroelectric phase transition using the effective- 
Hamiltonian approach of Ref. 



11. In order to obtain 



well-behaved smooth dielectric spectra especially closer 
to To in the low THz regime, MD time-step was cho- 
sen to be 0.5 fs and the simulations were conducted for 
a length of 13.2 ns in order to obtain 10000 samples for 
averaging the autocorrelation functions [t| . Figure 2 dis- 
plays the frequency dependence of the imaginary part of 
dielectric permittivity derived from MD simulations at 
270 K (C = 0.6 x 10 5 K and T = 245 K for simulated 
BST50). The presence of a lower frequency Debye-like 
CM in addition to the higher frequency SM is clearly 
seen there. 

To account for the simultaneously appearing SM and 
CM, we have used the coupled oscillator-relaxator model 
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FIG. 3: (Color online) Temperature dependence of the param- 
eters of Eq. 2 obtained from MD dielectric spectra of BST50. 
Filled symbols correspond to the unconstrained fit by Eq.[2] 
for temperatures below 400 K, full line in panel (b) stands 
for adjustment of the Arrhenius law (Eq.[9]) to 70 data below 
400 K. Open symbols were obtained for temperatures above 
Tcm from a constrained fit by Eq.[2]with 5 — 57 cm -1 and 70 
prescribed by Eq.[5]with U = 84meV and 700 = 1160 cm -1 . 



(following the notation of Ref.0): 



e(w) 



Sl§ 



7D<5 2 



(2) 



Here Sis, ^S, are the plasma frequency, oscillator fre- 
quency and damping constant of the SM, respectively; 
7d is the bare relaxation frequency of the CM, 6 is the 
coupling coefficient, and £co = 1 is the background per- 
mittivity [14|. Model parameters of the unconstrained 
fitting of Eq.[5]to the MD simulation data at various tem- 
peratures are shown as full symbols in Fig. 3. For tem- 
peratures above ~ 400 K, the CM is not visible in the 
simulated spectrum and the simple DHO formula (Eq.[2] 
with 5 = 0) provides a satisfactory fit (see Fig. 2). 

What can we learn from the temperature trends re- 
vealed in Fig. 3? First, let us note that Eq.[5]can be con- 
sidered as an asymptotic case of the dielectric response 
in the model of two coupled DHO's, one representing the 
normal phonon SM, and the other representing the ex- 
traneous excitation D [2l|: 



where the quantities S7d , , Td related to the CM have 
analogous meaning as Sis, wsj Ts- The coupling term A 
can be chosen positive. Indeed, in the Debye-relaxator 
limit (uj <C Trj and u -C wrj) the uj 2 term can be neglected 
against cj^ — iujTu and the above matrix reduces to that 
of the coupled relaxator-oscillator model Q . In this case, 
four parameters SId, Wd, Td and A can be replaced by 
only three parameters - the relaxational frequency 7d, 
the coupling frequency S, and an auxiliary frequency p 



7d = wn/r D , 5 = a/oju, p = w D n s /n D 



(5) 



With these substitutions (Eq.[5]) and for ui -C Tq, uj <C 
wd, and Sis 3> SId or p ^> 5, the expression of Eq.[3] 
reduces to the formula of Eq.[51 

Previously, several authors [1, 0, [B| treated the prob- 
lem of paraelectric CM in perovskites and derived a 
frequency dependence of the dielectric permittivity in 
the form equivalent to eqs.[2] or GHS Common assump- 
tion in all these approaches 0, [j, |f| is that the Curie- 
Weiss behavior of the static permittivity results from 
the Cochran-like temperature dependence of the SM fre- 
quency Wg = A(T — T c ) with a "Cochran" temperature 
T c < T ~ T c . This is indeed well obeyed here (with 
A = 50 cm^K" 1 and T c = 165 K for BST50, see Fig. 3). 
However, the implications for the temperature depen- 
dence of 8 and 7d are quite different among these models. 
In the spirit of Ref.0, one could assume that r5 and 7d 
either do not change with temperature since the bare D- 
mode is not contributing to the ordering mechanism, or 
perhaps, that (model A) 



and 



5 = const. 



7D ~ 1/T, 



(6) 



(7) 



when assuming in Eq. |(5j) that the D-mode damping in- 
creases with the temperature like Trj ~ T. Bare pseu- 
dospins in the model of Ref.0 are considered as freely 
fluctuating bistable dipoles forming an ideal paraelectric 
gas with Arrhenius relaxation law, which leads to (model 
B) 



a /T 



and 



7d = 7oo exp 



U 
kf 



(8) 



(9) 



(Uj) = £oo + [ Q S ] G(uj) 



with 

G(c)- 1 



— uj 2 — iuTu 
A* 



n s 



A 

,2 



ujg — uj — iujTs 



(3) 



, (4) 



where U is a potential barrier and 7oo is an attempt fre- 
quency. In a similar model applied to BaTiOs in Ref.[B|, 
one assumes that bare pseudospin subsystem would or- 
der at a finite temperature Tl < Tc, with Tl ~ Tc and 
the resulting formulas for 6 and 7d read (model C) 



S 2 = a /{T - Tl) 



(10) 



4 



and 



7d = 7oo exp 



U 
If 



T-T L 
T 



(11) 



Comparison of the anticipated temperature trends 
with the outcome of our MD simulations shown in Fig. 3 
indicates that none of the above three models is fully 
satisfactory. As a matter of fact, the MD results suggest 
that the coupling term 5 is constant (or it perhaps slightly 
increases with the increasing temperature), which is com- 
patible with the model A only, while the increasing 7^ 
rather supports the thermally activated flipping, assumed 
in models B or C. Tentatively, we have fitted the tem- 
perature dependence of 7d to the Arrhenius law of Eq.0 
in the 260-400 K temperature range. It gives reasonable 
values of U = 84meV and 70c = 1160 cm -1 . 

In fact, it turns out that the BST50 dielectric spectra 
can be well fitted at all temperatures (i.e., even above 
TcmO when keeping S = 57 cm -1 and 7d described by 
the aforementioned Arrhenius law (fitted parameters are 
shown by open symbols in Fig. 3). In other words, the 
coupled relaxator-oscillator model (Eq.[5]) can be used 
in the whole temperature range above Tc, without any 
particular discontinuity in the temperature course of the 
model parameters near the temperature Tcm- 

So why does the CM seem to "disappear" from the 
BST50 spectrum at Tom? In order to understand this 
behavior, it is convenient to consider the formula of Eq.[2] 
as the response of a DHO with a frequency-dependent 
self-energy term 



n(w) 



7 D 



r s 



<5 2 7d 
7d + w 



(12) 



From this expression, it is obvious that the simple 
DHO formula is recovered when 7d S> 6 2 /Ts, since 
in this case the imaginary part of the self-energy be- 
comes simply 7m[II] = — loTs- Therefore, it is reason- 
able to define the temperature Tcm as the temperature 
at which S 2 as 7dTs- Since the temperature dependence 
of 70 is steep, the temperature-driven crossover between 
7d -C 5 2 /T s and 7 D > S 2 /T s is quite sharp. With 
the data of Fig. 3, this crossover temperature is about 
Tcm ~ 360 — 370 K. Moreover, in case of Arrhenius law 
for 7d, one obtains 



T 



U 



cm 



fcln(7oor s /<5 2 ) 



(13) 



showing that Tcm scales with the activation energy U. 

This kind of crossover temperature Tcm likely bears 
some resemblance with the so-called Burns temperature 
and T* temperature, that have been observed in some 
complex perovskite materials HI- lil > since they all may 



be associated with change of dynamics. 

In summary, we provide interpretation of experiments 
and MD simulations suggesting that permittivity of BST 



system in the paraelectric phase shows a crossover be- 
tween a high-temperature spectrum with a simple SM 
and a lower-temperature spectrum with a more com- 
plex shape involving an additional CM that is coupled 
to the SM. Analysis of MD simulations allowed us to se- 
lect convenient formula for the temperature evolution of 
the model parameters, leading to a deviation with respect 
to previous theoretical predictions when CM and SM are 
coupled. As a matter of fact, these MD results indicate a 
thermally activated dynamics in a postulated pseudospin 
subsystem and their almost temperature-independent 
pseudospin-phonon coupling coefficient. Furthermore 
and interestingly, the dielectric spectra of BST50 from 
MD simulations could also be analyzed using this origi- 
nal coupled CM-SM model at any temperature above To 
(i.e., even in the high-temperature regime). As a result, 
the appearance of the CM in the dielectric spectra is 
understood as a crossover between fast and slow pseu- 
dospin dynamics. Arrhenius-type temperature depen- 
dence of bare pseudospin dynamics allows to understand 
the abrupt change of the spectrum at this crossover, 
which thus appears almost as a phase transition. 
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